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ABSTRACT 

Deep pencil beam surveys (< ldeg 2 ) are of fundamental importance for studying the high-redshift universe. 
However, inferences about galaxy population properties (e.g. the abundance of objects) are in practice limited 
by 'cosmic variance'. This is the uncertainty in observational estimates of the number density of galaxies 
arising from the underlying large-scale density fluctuations. This source of uncertainty can be significant, 
especially for surveys which cover only small areas and for massive high-redshift galaxies. Cosmic variance 
for a given galaxy population can be determined using predictions from cold dark matter theory and the galaxy 
bias. In this paper we provide tools for experiment design and interpretation. For a given survey geometry we 
present the cosmic variance of dark matter as a function of mean redshift z and redshift bin size Az. Using a 
halo occupation model to predict galaxy clustering, we derive the galaxy bias as a function of mean redshift for 
galaxy samples of a given stellar mass range. In the linear regime, the cosmic variance of these galaxy samples 
is the product of the galaxy bias and the dark matter cosmic variance. We present a simple recipe using a fitting 
function to compute cosmic variance as a function of the angular dimensions of the field, z, Az and stellar mass 
m*. We also provide tabulated values and a software tool. We find that for GOODS at z = 2 and with Az = 0.5 
the relative cosmic variance of galaxies with ra* > 10 u Mq is ~ 38%, while it is ~ 27% for GEMS and ~ 12% 
for COSMOS. For galaxies of m* ~ 1O 1O M the relative cosmic variance is ~ 19% for GOODS, ~ 13% for 
GEMS and ~ 6% for COSMOS. This implies that cosmic variance is a significant source of uncertainty at z = 2 
for small fields and massive galaxies, while for larger fields and intermediate mass galaxies cosmic variance is 
less serious. 

Subject headings: cosmology: theory — galaxies: high-redshift — galaxies: statistics — galaxies: stellar 
content — large-scale structure of universe 



1. INTRODUCTION 

Over the last decade, deep 'pencil beam' surveys of high- 
redshift galaxies have been the observational basis for study- 
ing the processes that drive galaxy formation and evolu- 
tion. For a given amount of observing time, such surveys 
have to trade off imaging area a nd imaging dept h: among 
HST imaging sur veys, GEMS dRix etalj 120041), AE GIS 
dDavis etalJ 120071) and COSMOS IScoville et all (20071) are 
examples that aim at comparatively wide fields to observe 
a larg e sample of galaxies. Alternatively, the Hubble Deep 
Field (IWilliams et al.|[T996l HDF) and the Ultra Deep Field 
dBeckwifh et alJl2006l UDF) are examples of extremely deep 
surveys in small areas, aimed at detecting faint galaxies 
(which can be either of low mass or at high redshift). 

One of the most fundamental properties of galaxy sub- 
populations at any epoch is their number density. However, 
observational estimates of galaxy number densities in finite 
volumes are subject to uncertainty due to cosmic variance, 
arising from underlying large-scale density fluctuations and 
leading to uncertainties in excess of naive Poisson errors. 
Note that this source of uncertainty is referred to as 'sample 
variance' in other branches of cosmology. For sampling vol- 
umes much larger than the typical clustering scale of the ob- 
served objects, cosmic variance is not significant. However, 
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many important existing surveys have a sampling volume that 
is small enough that cosmic variance may dominate the un- 
certainties. This may be particularly true at high redshift, 
where galaxies are expec ted to be much more strongly clus- 
tered than dark mat ter (Kauff mann et al.l 1 1 999fc [Baugh et al. 
11991 ICoil et "alHIOOlt Moster et al.H2009l) . Still, many pub- 
lished quantities which are based on number density (e.g. lu- 
minosity functions, stellar mass functions, etc.) are quoted 
with error budgets tha t do not properly accoun t for cosmic 
variance. As shown bv lTrenti~& Stiavell| d2008l) . the normal- 
ization and the slope of high-redshift luminosity functions can 
be affected by cosmic variance errors. 

Sound estimates of cosmic variance are helpful both for an- 
alyzing existing observational data and for designing future 
surveys. For a d etailed and general t r eatment of the cosmic 
err or, we refer tolSzapudi et alj dl999l) . IColombi et all d2000l) 
an d SzapudLeialJaOOO). Our goal here is to provide a simple 
recipe to derive cosmic variance for galaxy samples selected 
by their stellar mass m*, for a field of specified angular di- 
mensions a\ and «2. at a mean redshift z and for a redshift 
bin size of Az. Several previous works have provided esti- 
mates of co smic variance for speci fic sets of assumptions. For 
example, New man & Davis! (|2002) presented estimates of the 
cosmic variance as a function of field area and axis ratio for 
the redshift range 0.7 < z < 1.5, relevant to the DEEP2 red- 
shift survey, but did not account for galaxy bias (the cluster- 
ing amplitude of galaxie s relative to dark matter; see Eqn[T3b. 
Some rville et alj d2004l) provided predictions that could be 
used to estimate cosmic variance as a function of mean red- 
shift and survey volume, using the number density of the pop- 
ulation to estimate the bia s , assu ming one galaxy per halo. 
However, iSomerville et alj d2004l) used the approximation of 
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spherical (rather than re ctangular or pencil beam) volumes. 
iTrenti & Stiavellil ([2008) presented estimates of cosmic vari- 
ance for pencil beam geometries, again using the number den- 
sity of the population to estimate the host halo mass and hence 
the bias. 

Many galaxy samples in observational studies, however, are 
selected by their stellar mass, and stellar mass is one of the 
fundamental properties of galaxies. Galaxy clustering and 
bias are strong functions of stellar mass, making it extremely 
important to account for the mass dependence in cosmic vari- 
ance estimates. Moreover, one does not necessarily know the 
number density of a population of interest a priori. We there- 
fore take a different approach here, providing estimates of 
cosmic variance explicitly as a function of stellar mass and 
redshift. To do this, we make u se of the results previously 
presented in iMoster etaf] (120091 M09). We used a Halo Oc- 
cupation Distribution (HOD) model to empirically establish 
the relationship between stellar mass and dark matter halo (or 
sub-halo) mass at different redshifts, and then used a dissi- 
pationless N-body simulation to compute the galaxy bias as 
a function of stellar mass and redshift. We presented fitting 
functions for this quantity in M09. We combine these with 
estimates of the cosmic varia nce for the dark matter f or rect- 
angular cells as described in Newman & Davis (2002). 

The paper is organized as follows: in Section [2] we de- 
scribe how to define and how to compute cosmic variance for 
a galaxy population, based on the underlying cold dark matter 
(CDM) theory, a model for galaxy bias and for a combination 
of separate survey fields. In Section [3] we spell out a recipe 
to compute cosmic variance for a pencil beam geometry in 
four steps: (1) select survey geometry, (2) choose mean red- 
shift and redshift bin size, (3) determine stellar mass interval 
and (4) compute cosmic variance. We also apply this recipe 
to four existing surveys at a redshift of z = 2 and z = 3.5. In 
Section H] we compare cosmic variance for different galaxy 
samples at different redshifts and different survey geometries. 
Finally, we summarize our methods and conclusions in Sec- 
tion 

Throughout, we assume cosmo logical parameters c onsis- 
tent with results from WMAP-3 fepergel e tlri] 120071) for a 
flat ACDM cosmological model: matter density fi m = 0.26, 
cosmological constant £7a = 0.74, Hubble parameter Hq = 
72 km s -1 Mpc -1 , fluctuation amplitude erg = 0.77 and primor- 
dial power -spectrum n s = 0.95. All stellar masses are cali- 
brated to a Kroupa (2001) initial mass function. 

2. METHOD 

The mean (N) and the variance (N 2 )-(N) 2 are given by 
the first and second moments of the probability distribution 
Pn(V), which describes the probability of counting N objects 
within a volume V. The relative cosmic variance is defined as 



be written as 



{N 2 } _ {N) 2_ {N) 
W 2 



(1) 



The last term in the numerator represents the correction for 
Poisson shot noise. The second moment of the object counts 
is 



(AO! 

V 2 



dV,dVb£(|r a -rfe| 



(2) 



where £ is the two-point correlati on function of the sample 
and V is the sample volume (see Peebles 1980, pg. 234). 
Combining this with Equation (Q}, the cosmic variance can 
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(3) 



There are two approaches to evaluating this equation. The 
first is applicable to populations with a known correlation 
function, when the integral can be solved either analytically 
or numerically as a function of cell radius or cell volume. 
However, in practice, there are a number of difficulties with 
this approach: while on small scales in the nonlinear regime 
(r < 10- 15 Mpc) the correlation function can typically be ap- 
proximated by a power law, it deviates on larger scales (in the 
linear regime), from this form. Additionally, in many cases of 
interest the correlation function is not known a priori. 

The second approach to calculating the cosmic variance can 
be used if the correlation function cannot be approximated by 
a power law or if it is unknown. It makes use of the galaxy 
bias, which can be predicted for a given galaxy population 
using halo occupation models. 

We can substitute the galaxy correlation function £ gg in 
Equation (0) for the correlation function for dark matter 



£ gg (r,m,z) = b (m*,z) & m (r,z) , 



(4) 



where b(m*,z) is the galaxy bias, which depends both on red- 
shift z, and the stellar mass of the galaxies m*. In general, 
the bias is also a function of r, however, in the linear regime, 
we assume that it is independent of scale. This assumption 
will be made throughout the paper. Employing the bias in this 
manner yields 

tr 2 (m„,z)=^2 J dV R dV b b 2 (m*,z) &/ m (r ab ,z) 



--b 2 (m*,z) ^ / AYmYu ;,„„(>>.;.) 
--b 2 (m^,z) oj m (z) 



(5) 



The integral leading to <jj m can then be solved using linear 
theory in CDM for any volume V with a given geometry. The 
cosmic variance for galaxies can thus be determined by mul- 
tiplying the dark matter cosmic variance at a given redshift 
with the galaxy bias for a given stellar mass at that redshift. 
For a pencil beam geometry the volume element V is set by 
the angular size of the survey and the redshift bin size. 

We now pursue this second approach, using the galaxy bias 
as a function of re dshift and stellar mass as computed by 
Moster et al. (2009). In this paper the authors used a halo 
occupation model to obtain a parameterized stellar-to-halo 
mass relation by populating halos and subhalos in an iV-body 
simulation with galaxies and requiring that the observed stel- 
lar mass function be reproduced. Using the halo positions 
obtained from the simulation, they find that predictions for 
the galaxy correlation function are in excellent agreement 
with observed clustering properties at low redshift. The de- 
rived stellar-to-halo mass relation is finally used to predict the 
stellar-mass-dependent galaxy correlation function and bias at 
high redshift. 

We calculate the cosmic variance for dark matter in each 
redshift bin by solving the integral in Eqn (O for dark 
matter using the code Q UICKCVP] which is described in 
iNewman & Davisl d2002). The code computes the dark mat- 
ter correlation function using the transfer function given in 

1 QUICKCV is available at http://astro.berkeley.edu/~jnewman/research 
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iBardeen et al.l d 19861) and a window function which is 1 in- 
side the volume V and elsewhere. The volume elements are 
determined by the angular dimensions of the survey field and 
the redshift bin. Finally, the cosmic variance for galaxies is 
computed using Equation ©. 

Once cosmic variance has been determined for a number of 
individual survey fields that are widely separated on the sky, 
compared to the correlation length, as is often the case (e.g. 
for the various HST survey fields), one can easily derive the 
cosmic variance for the combined sample. For the combina- 
tion of fields i the total volume is 



V t0 , = E, V t , 

implying a cosmic variance for the total sample of 

<& = W i dVadV b a|ra-r b |) 
not Jv 

= 7^2 / d(E ! V0ad(E;V 7 ) b £(|r, a -i>|) 

[X-i Vi) JV 

/ dV^dV^Clria-r^l) 
i Vi) Jv 

^—s E, / dy ;a dy ;b a|r ;a -r ;b |) 

A Vi) Jv 



(6) 



(EfVi) 

(eT 



(7) 



(Si Vd l 

where we have used 

a|r, a -i>|) = 0for*yj (8) 

since the galaxies in separated fields are uncorrected. For a 
pencil beam geometry (angular extend <C 1 radian) this can 
be written as 

2 S,-(aia 2 )?cr, 2 , m 

°tot = 2 W 

[S; (aia 2 )i] 

where a\ and a 2 are the angular dimensions of each field. 
Note that this volume-weighted average is identical to an 
inverse-variance-weighted average only in the range where 
the cosmic variance scales with volume. 

3. RECIPE FOR COSMIC VARIANCE 

In order to compute the cosmic variance for galaxies in a 
given volume in the universe one has to fix five parameters. 
The size of the volume is determined by two dimensions tan- 
gential to the line of sight and one dimension parallel to the 
line of sight. The tangential dimensions we use are the angu- 
lar size of the field a\ and a 2 , while the vertical dimension is 
set by the size of the redshift bin under consideration Az. The 
fourth parameter is the mean redshift of the volume z, since 
the cosmic variance depends on the clustering strength which 
is a function of redshift. Moreover z determines the length 
scale of the volume bins by setting the conversion factor from 
angular dimensions and Az into comoving Mpc. In order to 
determine the cosmic variance for galaxies one finally has to 
fix the stellar mass of the sample under consideration, as 
the clustering strength and the bias depend strongly on stellar 
mass. Thus the five parameters which have to be set are ol\, 
a>2, z, Az and ra*. 

3.1. Survey Geometry 

The first step to calculate the cosmic variance is to fix the 
size of the volume element under consideration. For now 



Table 1 

Summary of the survey geometries 



Survey 






a 2 " 


area fle id b 


total area b 


UDF 


1 


3.3 


3.3 


11 


11 


GOODS 


2 


10 


16 


160 


320 


GEMS 


1 


28 


28 


784 


784 


EGS 


1 


10 


70 


700 


700 


COSMOS 


1 


84 


84 


7056 


7056 



L Angular dimension in arcmin 
' Areas in arcmin 2 




UDF A A 

GOODS *- X 

GEMS O O 

EGS + + 

COSMOS X X 



A a 



Aa-a-a-A 



12 3 

z 

Figure 1. Cosmic variance as a function of mean redshift with a reference b in size of 
Az re f = 0.2. The different sym b ols correspond to different fields while the lines are fits 
to the symttols using Equation 1101 . 

we set the size of the redshift bin to a reference value of 
Az = 0.2. This is a convenient bin size for low redshift while 
for higher redshift it is rather small. However, this choice can 
be changed at any redshif t (fo r values Az > 0.05) using the 
recipe presented in section [3~2l For the tangential dimensions, 
we choose five popu lar fields: the Hubble Ultra Deep Field 
dBeckwifh et alj|2006l UDF X one field of the Grea t Observa- 
tories Origins Deep Survey dGiavalisco et al.ll2004 GOODS), 
the Extended Chandra Deep Field South (E-CDFS) used in 
the Galaxy Evo lution From Morphology And SEDs survey 
dRix etalJl2004l GEMS), the Extended Groth Strip (EGS) 
used in t he All-wavelength Extended Groth strip International 
Survey dDavis et all 120071 AEGIS) and the field used in the 
Cosmic Evolution Survey dScoville etalj 120071 COSMOS). 
The angular dimensions of the surveys and the resulting ob- 
servation areas are summarized in Table Q] We note that due 
to the small scales for the UDF numerical uncertainties in the 
integration of the correlation function are likely to be high. As 
a result the values for cosmic variance in the UDF computed 
here, especially for z < 0.5 (where size of the volume ele- 
ments is very small) have considerable uncertainties. At these 
scales, the assumption of linear biasing is also not accurate. 
As a consequence the cosmic variance for the UDF presented 
here is to be treated with caution and is mainly included to 
demonstrate the large cosmic variance in small fields. 
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Table 2 

Cosmic variance for different surveys with common redshift bins. The first 
and second columns show the mean redshift and the redshift bin size, the 

third column gives the root comic variance for dark matter while the last six 
columns give the root cosmic variance for galaxies in the indicated stellar 
mass intervals. 
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<Tgg for log(m*/M )±O.25 









8.75 


9.25 


9. / J 


1U.ZJ 


in n ^ 

1 U. 1 J 
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UDF 3.3' 


x 3.3 








0.1 


0.2 


0.316 


0.349 


0.358 


0.380 


0.409 


0.430 


0.526 


0.3 


0.2 


0.237 


0.272 


0.288 


0.295 


0.318 


0.339 


0.428 


0.5 


0.2 


0.218 


0.262 


0.273 


0.283 


0.307 


0.332 


0.435 


0.7 


0.2 


0.204 


0.259 


0.271 


0.280 


0.305 


0.335 


0.456 


0.9 


0.2 


0.192 


0.259 


0.274 


0.282 


0.308 


0.346 


0.488 


1.1 


0.2 


0.182 


0.264 


0.281 


0.289 


0.317 


0.365 


0.531 


1.5 


0.2 


0.167 


0.282 


0.305 


0.319 


0.352 


0.426 


0.648 


1 


u.z 


ft 1« 
U. 1 JO 


0.311 


0.341 


U. JOV 




U. J 1 a 


U.oU / 


2.5 


0.2 


0.143 


0.374 


0.416 


0.482 


0.529 


0.718 


1.129 


3.5 


0.2 


0.128 


0.523 


0.593 


0.781 


0.845 


1.236 


1.898 








GOODS 1C 


i' x 16' 








0.1 


0.2 


0.244 


0.270 


0.277 


0.294 


0.317 


0.333 


0.407 


0.3 


0.2 


0.205 


0.236 


0.243 


0.255 


0.276 


0.293 


0.371 


0.5 


0.2 


0.179 


0.215 


0.224 


0.232 


0.252 


0.272 


0.357 


0.7 


0.2 


0.160 


0.203 


0.213 


0.220 


0.239 


0.263 


0.358 


0.9 


0.2 


0.146 


0.197 


0.208 


0.214 


0.234 


0.263 


0.371 


1.1 


0.2 


0.135 


0.195 


0.207 


0.214 


0.235 


0.270 


0.392 


1.5 


0.2 


0.119 


0.200 


0.217 


0.227 


0.250 


0.303 


0.460 


1.9 


0.2 


0.107 


0.215 


0.235 


0.255 


0.281 


0.357 


0.557 


2.5 


0.2 


0.096 


0.250 


0.278 


0.322 


0.354 


0.480 


0.754 


3.5 


0.2 


0.083 


0.336 


0.382 


0.503 


0.544 


0.795 


1.221 










GEMS 28' 


X28' 









0.1 


0.2 


0.226 


0.250 


0.256 


0.272 


0.293 


0.308 


0.377 


0.3 


0.2 


0.174 


0.200 


0.206 


0.216 


0.234 


0.249 


0.315 


0.5 


0.2 


0.144 


0.173 


0.180 


0.187 


0.203 


0.219 


0.288 


0.7 


0.2 


0.124 


0.158 


0.165 


0.171 


0.186 


0.204 


0.279 


0.9 


0.2 


0.110 


0.149 


0.157 


0.162 


0.177 


0.199 


0.280 


1.1 


0.2 


0.100 


0.144 


0.154 


0.158 


0.174 


0.200 


0.291 


1.5 


0.2 


0.085 


0.143 


0.155 


0.162 


0.179 


0.217 


0.329 


1.9 


0.2 


0.075 


0.150 


0.164 


0.178 


0.196 


0.250 


0.389 


2.5 


0.2 


0.065 


0.170 


0.189 


0.219 


0.240 


0.326 


0.512 


3.5 


0.2 


0.054 


0.221 


0.251 


0.331 


0.358 


0.523 


0.804 










EGS 10' 


x 70' 








0.1 


0.2 


0.215 


0.238 


0.244 


0.259 


0.279 


0.293 


0.358 


0.3 


0.2 


0.157 


0.180 


0.186 


0.195 


0.211 


0.224 


0.284 


0.5 


0.2 


0.128 


0.154 


0.160 


0.167 


0.181 


0.195 


0.256 


0.7 


0.2 


0.120 


0.140 


0.146 


0.151 


0.164 


0.181 


0.246 


0.9 


0.2 


0.097 


0.131 


0.139 


0.143 


0.156 


0.175 


0.247 


1.1 


0.2 


0.088 


0.127 


0.135 


0.139 


0.153 


0.176 


0.256 


1.5 


0.2 


0.075 


0.126 


0.136 


0.143 


0.157 


0.191 


0.290 


1.9 


0.2 


0.066 


0.132 


0.145 


0.157 


0.173 


0.220 


0.343 


2.5 


0.2 


0.057 


0.150 


0.167 


0.194 


0.213 


0.288 


0.453 


3.5 


0.2 


0.048 


0.197 


0.224 


0.295 


0.319 


0.466 


0.716 



COSMOS 84' X 84' 



0.1 


0.2 


0.181 


0.200 


0.205 


0.218 


0.235 


0.247 


0.302 


0.3 


0.2 


0.117 


0.134 


0.138 


0.145 


0.157 


0.167 


0.211 


0.5 


0.2 


0.088 


0.105 


0.109 


0.114 


0.123 


0.133 


0.175 


0.7 


0.2 


0.070 


0.089 


0.094 


0.097 


0.105 


0.116 


0.158 


0.9 


0.2 


0.059 


0.080 


0.084 


0.087 


0.095 


0.107 


0.151 


1.1 


0.2 


0.051 


0.074 


0.079 


0.083 


0.090 


0.103 


0.150 


1.5 


0.2 


0.041 


0.070 


0.075 


0.079 


0.087 


0.105 


0.160 


1.9 


0.2 


0.035 


0.070 


0.077 


0.083 


0.092 


0.117 


0.182 


2.5 


0.2 


0.029 


0.076 


0.085 


0.098 


0.108 


0.146 


0.230 


3.5 


0.2 


0.023 


0.095 


0.108 


0.142 


0.154 


0.225 


0.346 



Table 3 

Fitting parameters for different surveys 



Survey 



P 



UDF 
GOODS 
GEMS 
EGS 

COSMOS 



0.251 
0.261 
0.161 
0.128 
0.069 



0.364 
0.854 
0.520 
0.383 
0.234 



0.358 
0.684 
0.729 
0.673 
0.834 



3.2. Redshift bins 

The next step is to calculate the cosmic variance for dark 
matter for a set of mean redshifts z. We do this using the 
ang ular dimensions of the five surveys presented in section 
3. H and the reference redshift bin size of Az = 0.2. The cosmic 
variance is then computed for a set of mean redshifts. The 
resulting values of the root cosmic variance for dark matter 
are given in the third column of Table [2] and are plotted in 
Figure[T](symbols). To make it possible to obtain the cosmic 
variance for mean redshifts other than these tabulated values, 
we introduce a fitting function: 



<Jdm(z, Az = 0.2) : 



Ob 



(10) 



The parameters a a , at, and j3 depend on the angular dimen- 
sions of the field and are given in Table [3] The lines in Figure 
Q]show the fitting function for the five surveys and agree very 
well with the computed values for the redshift bins. Using 
Equation (fTOb for any field we are thus able to calculate the 
dark matter cosmic variance as a function of mean redshift. 

Up until now we have assumed a reference redshift bin size 
of Az re f = 0.2 in order to determine the size of the volume 
element under consideration. However, especially at higher 
redshift, larger bin sizes are desirable. We thus investigate 
how the cosmic variance depends on Az. For this we plot 
(Tdm(Az)/o'dm(Az r ef) as a function of Az/Az re / for different 
fields and mean redshifts in Figure [2] We find that indepen- 
dent of angular dimensions and mean redshift the cosmic vari- 
ance has the same dependance on the redshift bin size: 



,(Az) / Az 



f7 d m(Az re f) \Az r 



ef 



(ID 



This means that for a reference redshift bin size of Az re f = 0.2 
as assumed in Equation ( TTOt the root cosmic variance can be 
calculated for a different Az as: 



6T d m(Az,z) = o- dm (Az = 0.2, z) 




(12) 



There is a simple picture that explains the dependence on Az. 
If we divide a redshift bin of size Az into TV redshift bins 
with size Az, then Az = NAz,. Assuming that all these red- 
shift bins are uncorrelated, we can invoke Equation (O which 
yields er(Az) = cr(Azi)/ \/N- Using TV = Az/ Azi we finally get 
er(Az) = cr(Azi) x yj Az,/ Az. It is interesting that the sim- 
plified assumption of uncorrelated volume elements yields a 
relation that fits so well to the results calculated for the total 
redshift bin, as seen in Figure [2] However, we note that this 
approximation is only valid for redshift bin sizes Az > 0.05. 
For smaller bin sizes w e refer to the mor e accurate model pre- 
sented in section 3.5 of [Newman] ((2008). 
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Table 4 

Galaxy bias fit parameters 



logm ? 


bo 


bi 


b 2 


8.5 


- 9.0 


0.062 ±0.017 


2.59 ±0.18 


1.025 ±0.062 


9.0 


- 9.5 


0.074 ± 0.008 


2.58 ± 0.26 


1.039 ±0.028 


9.5 


- 10.0 


0.042 ± 0.003 


3.17 ± 0.05 


1.147 ±0.021 


10.0 


- 10.5 


0.053 ± 0.014 


3.07 ±0.17 


1.225 ±0.077 


10.5 


- 11.0 


0.069 ± 0.014 


3.19 ±0.13 


1.269 ±0.087 


11.0 


- 11.5 


0.173 ±0.035 


2.89 ± 0.20 


1.438 ±0.061 


> 


8.5 


0.063 ± 0.008 


2.62 ± 0.08 


1.104 ±0.028 


> 


9.0 


0.085 ± 0.008 


2.50 ± 0.06 


1.098 ±0.064 


> 


9.5 


0.058 ± 0.005 


2.96 ± 0.06 


1.192 ±0.021 


> 


10.0 


0.072 ±0.013 


2.90 ±0.13 


1.257 ±0.051 


> 


10.5 


0.093 ± 0.015 


3.02 ±0.11 


1.332 ±0.054 


> 


11.0 


0.185 ±0.032 


2.86 ± 0.25 


1.448 ±0.098 



Note. — All quoted masses are in units of Mq 

3.3. Stellar mass dependence 

The last step to derive the cosmic variance for a given sam- 
ple of galaxies is to apply the galaxy bias. As we have shown 
in Equation (O, in the linear-biasing limit, the cosmic vari- 
ance for a sample of galaxies with a given stellar mass is the 
product of the squared bias and the dark matter cosmic vari- 
ance. The stellar mass dep endent galaxy bias has been de- 
rived in Moster et al. (2009). The authors present the redshift 
dependence of the bias using parameterized functions of the 
form: 

b(m*,z) = b (z+l) bl +b 2 (13) 

with the parameters bo, b\, and b 2 which are given in Table 
2]for six stellar mass bins and six stellar mass thresholds. In 
order to compute the root cosmic variance for galaxies of that 
mass, we have to multiply the bias and the root cosmic vari- 
ance of dark matter for the same redshift. This has been done 
for the five surveys presented in Table [2] the cosmic variance 
for galaxies of mass m* is given in columns 4-9 for the refer- 
ence redshift bin size of Az = 0.2. Cosmic variance for mas- 
sive galaxies is always larger than that of low mass galaxies, 
since the galaxy bias increases with increasing stellar mass. 




3.4. Cookbook for Cosmic Variance 

Finally we can summarize our recipe to derive the cosmic 
variance for a particular survey: 

Step 1): Choose the survey (field) in Table[T]that is closest to 
the field you are using. 

Step 2): Choose a mean redshift z and a redshift bin size Az. 

Step 3): Choose the stellar mass range of your galaxy sample 
from Table g] 

Step 4a): Calculate the dark matter root cosmic variance 
ad m (z, Az = 0.2) using Equation ( fT0T > and the parame- 
ters for your survey as given in Table [3] 

Step 4b): Calculate the galaxy bias b(m*,z) using Equation 
(TPTt and the parameters for your stellar mass bin or 
threshold as given in Table |4] 

Step 4c): Compute the root cosmic variance for galaxies: 

o-gg(m,,z, Az) = fr(m»,z) a dm (z, Az = 0.2) 
Step 4d): Combine widely separated fields using Equation 

This recipe allows for easy computation of cosmic variance 
for the representative surveys we have presented. To enable 
readers to compute cosmic variance for different survey areas 
or geometries, we provide a software tool as an online sup- 
plementQ. For specified angular dimensions and redshift bins 
the code computes the cosmic variance for dark matter and 
galaxies in six stellar mass bins. 

3.5. Examples 

Having derived the recipe to determine the cosmic vari- 
ance, we can apply it to the five example surveys. We first 
give an overview of the surveys. The UDF is an exposure 
of a square field in the southern sky covering a total area of 
1 1 arcmin 2 with the Advanced Camera for Surveys (ACS) on 
the Hubble Space Telescope (HST). To date it is the deep- 
est image of the universe with uniform limiting magnitudes 
m ~ 29 for point sources and contains at least 10,000 ob- 
jects. However, of the five surveys studied here, the UDF 
is the smallest. The GOODS project covers two fields, the 
Hubble Deep Field North (HDFN) and the Chandra Deep 
Field South (CDFS) which both have similar dimensions of 
« 10 x 16 arcmin 2 and are so widely separated as to be uncor- 
rected. Since for widely separated fields the cosmic variance 
goes as 1 /Vfieid, the total variance for both fields decreases by 
a factor of 2. The GEMS survey has imaged the E-CDFS, a 
square area of 784 arcmin 2 , with the ACS on the HST. The E- 
CDFS is centered on the CDFS and contains roughly 10,000 
galaxies down to a depth of 24th magnitude in the R-band. 
AEGIS is a multi-wavelength, deep, wide-field, photometric 
and spectroscopic survey in the Extended Groth Strip (EGS) 
area. While the total area is similar to that of GEMS, the field 
is a long strip with dimensions of w 10 x 70 arcmin 2 . The 
widest survey for which we compute cosmic variance is COS- 
MOS which covers an equatorial square field of 2 deg 2 with 

2 IDL code available at http://www.mpia.de/homes/moster/research 
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Figure 3. The root cosmic variance as a function of redshift with redshift bins of 
constant volume. The four panels are for dark matter and galaxies in three stellar mass 
bins. The lines are for different surveys. 

imaging by most of the major space-based telescopes and sev- 
eral large ground-based telescopes. The magnitude limit in 
the I-band is / — 27. 

If we apply this recipe to the UDF, GOODS, GEMS, 
AEGIS and COSMOS at z = 2 and Az = 0.5, the relative 
cosmic variance for massive galaxies with m* > 10 u Mq is 
~ 55% for the UDF, while it is - 38% for GOODS, - 27% 
for GEMS, - 23% for AEGIS and - 12% for COSMOS. 
For galaxies of intermediate mass (m* ~ 10 10 M Q ) the relative 
cosmic variance is ~ 27% for the UDF, ~ 19% for GOODS, 
- 13% for GEMS, - 11% for AEGIS and - 6% for COS- 
MOS. This implies that at z = 2 cosmic variance is a signifi- 
cant source of uncertainty for small fields and massive galax- 
ies. However, for larger fields and intermediate mass galaxies 
cosmic variance is less serious. 

At a higher redshift of z = 3 .5 and Az = 0.5 the relative cos- 
mic variance is higher, since the galaxy bias is much larger. 
For massive galaxies with m* > 1O U M , we find a relative 
cosmic variance of - 124% for the UDF, ~ 78% for GOODS, 
~ 51% for GEMS, - 45% for AEGIS and - 21% for COS- 
MOS. For intermediate mass galaxies with — 1O 1O M the 
relative cosmic variance is — 54% for the UDF, — 34% for 
GOODS, ~ 22% for GEMS, - 20% for AEGIS and - 9% 
for COSMOS. This shows that at z = 3.5 cosmic variance can 
be serious in small fields even for intermediate mass galaxies. 
For massive galaxies even the widest survey (COSMOS) has 
considerable uncertainties due to cosmic variance. 

4. SCALING OF COSMIC VARIANCE 

It is not very useful to compare the cosmic variance at dif- 
ferent redshifts for a fixed bin size, because the volumes are 
not the same. Since cosmic variance strongly depends on the 
volume of the sample under consideration, we choose redshift 
bins such that all redshift bins have the same volume. This al- 
lows us to directly compare the cosmic variance at different 
redshifts. The results of this analysis are given in Table[5]for 
the GOODS, GEMS and COSMOS surveys. 



Table 5 

Cosmic variance for different surveys with constant comoving volume. The 
first and second columns show the mean redshift and the redshift bin size, 
the third column gives the root comic variance for dark matter, and the last 
six columns give the root cosmic variance for galaxies in the indicated 
stellar mass intervals. 
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8.75 


9.25 


9.75 


10.25 


10.75 


11.25 








GOODS 10' 


x 16' 








0.00 


1.12 


0.081 


0.099 


0.103 


0.107 


0.116 


0.126 


0.168 


1.12 


1.58 


0.084 


0.133 


0.143 


0.148 


0.163 


0.194 


0.291 


1.58 


1.99 


0.078 


0.149 


0.163 


0. 174 


0. 192 


0.241 


0.374 


1.99 


2.39 


0.073 


0.165 


0.183 


0.204 


0.225 


0.295 


0.464 


2.39 


2.78 


0.069 


0.186 


0.208 


0.243 


0.267 


0.365 


0.573 


2.78 


3.17 


0.065 


0.209 


0.235 


0.290 


0.316 


0.446 


0.696 


3.17 


3.58 


0.060 


0.231 


0.262 


0.339 


0.368 


0.534 


0.823 


3.58 


4.00 


0.056 


0.260 


0.296 


0.404 


0.436 


0.647 


0.983 








GEMS 28' 


x 28' 








0.00 


1.12 


0.063 


0.077 


0.080 


0.083 


0.091 


0.098 


0.130 


1.12 


1.58 


0.060 


0.096 


0.103 


0.107 


0.118 


0.140 


0.210 


1.58 


1.99 


0.055 


0.105 


0.114 


0.122 


0.135 


0.169 


0.262 


1.99 


2.39 


0.050 


0.114 


0.126 


0.140 


0.155 


0.203 


0.319 


2.39 


2.78 


0.047 


0.126 


0.141 


0.165 


0.181 


0.247 


0.389 


L. to 


j.L 1 


U.U43 


0.140 


0.158 


0.194 


0.212 


0.299 


U.40D 


3.17 


3.58 


0.040 


0.153 


0.173 


0.225 


0.244 


0.354 


0.545 


3.58 


4.00 


0.037 


0.171 


0.194 


0.265 


0.286 


0.425 


0.645 








COSMOS 84 


' X 84' 








0.00 


1.12 


0.036 


0.044 


0.046 


0.048 


0.052 


0.057 


0.075 


1.12 


1.58 


0.030 


0.047 


0.051 


0.053 


0.058 


0.069 


0.104 


1.58 


1.99 


0.026 


0.049 


0.054 


0.058 


0.064 


0.080 


0.124 


1.99 


2.39 


0.023 


0.052 


0.057 


0.064 


0.070 


0.093 


0.145 


2.39 


2.78 


0.021 


0.056 


0.063 


0.074 


0.081 


0.110 


0.173 


2.78 


3.17 


0.019 


0.061 


0.069 


0.085 


0.093 


0.131 


0.204 


3.17 


3.58 


0.017 


0.066 


0.075 


0.097 


0.105 


0.153 


0.235 


3.58 


4.00 


0.016 


0.073 


0.083 


0.113 


0.122 


0.181 


0.275 



4.1. Cosmic variance for different redshifts 

We would expect the cosmic variance of dark matter to de- 
crease monotonically with increasing redshift, because it de- 
pends on the clustering which is lower for high redshift. How- 
ever, we see, that for the UDF adm increases, reaches a max- 
imum and then decreases with increasing redshift, while for 
the other surveys adm decreases with increasing redshift. The 
reason for this is that although the different redshift bins have 
the same comoving volume, their geometry is not the same. 
This is due to the pencil beam geometry for which the volume 
of the low redshift bins are long but narrow, while the higher 
redshift bins become shorter and wider. Thus with increasing 
redshift the ratio of the area and the depth of the redshift bins 
increases. In the lower redshift bins the galaxies are thus on 
average more separated than in the higher redshift bins. This 
causes the cosmic variance to increase with redshift. At high 
redshift and for very wide surveys the ratio decreases again 
with redshift as the line of sight distance in a redshift bin be- 
comes smaller than the transverse distance. 

Both effects, the decrease of adm with redshift due to lower 
clustering amplitudes and the growth of adm because of the 
geometry of the redshift bins, affect ad m (z)- The redshift at 
which ad„, is maximal depends on the area of the survey. For 
small areas (e.g. UDF) it is higher (z max —1.3) while for large 
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Figure 4. Effects of survey geometry on cosmic variance. The root cosmic variance 
of dark matter is plotted as a function of the ratio between the two observation angles 
a\ j ol2- The lines represent different redshift bins with a bin size of Az = 1 .0. 

areas (e.g. COSMOS) it is lower (Zmax < 0.2). The cosmic 
variance for the galaxy samples we consider, however, always 
increases with redshift. This is due to the large biases of these 
samples, which are always larger than unity. Since the bias in- 
creases strongly with redshift, the cosmic variance for galax- 
ies thus also increases with redshift. 

Figure [3] plots the cosmic variance as a function of red- 
shift with redshift bins of constant volume. The four different 
panels show cr v for dark matter and three stellar mass bins. 
The lines are for different surveys. We see that except for 
the UDF, the cosmic variance for dark matter decreases with 
redshift while for all surveys the cosmic variance for galaxies 
increases with redshift. 

4.2. Cosmic variance for different geometries 

As we have shown, cosmic variance depends strongly on 
the geometry of the survey. In section |4~T1 we investigated the 
geometric effect that arises because of the ratio between line- 
of-sight distance and the transverse distance in the redshift 
bins. Now we investigate the effects that arise due to different 
angular geometry on the sky. For this we assume a survey 
with a fixed total area of l.Odeg 2 . We vary the ratio between 
the two observation angles a\ joti from a\/a2 ~ to a\/ct2 = 
1 . This is done for several mean redshifts with a bin size of 
Az = 0.5. 

Figure|4]shows the geometry effects on the cosmic variance. 
It plots o~dm as a function of the ratio between the two observa- 
tion angles. The result agrees with our explanation in section 
14.11 for a low ratio the cosmic variance is very low. If one 
increases the ratio the cosmic variance quickly increases and 
reaches its maximum at a \jct2 = 1. Here the mean distance 
between the observed galaxies is the smallest, so that many 
galaxies are likely to be correlated, resulting in large cosmic 
variance. 

This effect can also be seen by comparing the cosmic vari- 
ance for GEMS and EGS at a fixed redshift bin. We find 
that although the area of the GEMS field is larger than that 
of the EGS, GEMS has a higher cosmic variance. However, 
we note that unless the ratio between the observation angles 
is smaller than ~ 20%, th is effect is small. In agreement with 
iNewman & Davisl ((2002), we find that the size of the survey 



area is much more important than the axis ratio. 

5. CONCLUSIONS 

Inferences about galaxy number densities and related quan- 
tities such as luminosity and stellar mass functions are subject 
to uncertainties due to cosmic variance. It can be a significant 
source of uncertainty, especially for surveys which cover only 
small areas and for massive high-redshift galaxies. We have 
derived a simple recipe to compute cosmic variance for five 
surveys as a function of mean redshift z, redshift bin size Az 
and the stellar mass of the galaxy population to». 

For this recipe we first calculated the dark matter cosmic 
variance as a function of mean redshift and for a reference 
bin size of Az re d = 0.2. This was done by integrating the dark 
matter correlation function obtained from CDM theory for a 
pencil beam geometry in a redshift bin. We provide fitting 
functions for five surveys (UDF, GOODS, GEMS, AEGIS, 
COSMOS). We then showed that the dependence of cosmic 
variance on bin size Az is nearly independent of mean redshift 
and survey geometry. This means that cosmic variance for a 
different bin size can be estimated from the value computed 
with the reference bin size using a simple co nversion fac- 
tor. W e have used the galaxy bias predictions by Most er et al.l 
(2009) (for six stellar mass bins and six thresholds) and the 
dark matter cosmic variance to compute the cosmic variance 
for galaxies. We also presented a formula to compute the 
overall cosmic variance for a volume-weighted average of 
multiple separated fields. 

Applying this recipe to GOODS, GEMS and COSMOS at 
z = 2 and Az = 0.5, the relative cosmic variance for galax- 
ies with to, > 10 n M© (to, - 10 10 M Q ) is ~ 38% (19%) for 
GOODS, while it is ~ 27% (13%) for GEMS and - 12% (6%) 
for COSMOS. Atz = 3.5 and Az = 0.5 we found a relative cos- 
mic variance for galaxies with to, > 10 n M Q (m, ~ 1O 1O M ) 
of - 78% (34%) for GOODS, - 51% (22%) for GEMS and 
-21% (9%) for COSMOS. 

We find that for all fields the cosmic variance is much more 
significant for massive galaxies, which is due to the strong de- 
pendence of the galaxy bias on stellar mass. Our results imply 
that cosmic variance is a significant source of uncertainty at 
Z = 2 for small fields and massive galaxies, while for larger 
fields and intermediate mass galaxies, cosmic variance is less 
serious. At z = 3.5 cosmic variance can become serious in 
small fields even for intermediate mass galaxies. For massive 
galaxies even the widest survey (COSMOS) has significant 
uncertainties due to cosmic variance. 

BPM thanks the Space Telescope Science Institute for hos- 
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